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Abstract 

The top-quark mass is measured in proton-proton collisions at ^Js = 7 TeV using a 
data sample corresponding to an integrated luminosity of 5.0 fb -1 collected by the 
CMS experiment at the LHC. The measurement is performed in the dilepton decay 
channel tt — > (£ + Veb) (£~Vgb), where i = e,]i. Candidate top-quark decays are se- 
lected by requiring two leptons, at least two jets, and imbalance in transverse mo- 
mentum. The mass is reconstructed with an analytical matrix weighting technique 
using distributions derived from simulated samples. Using a maximum-likelihood 
fit, the top-quark mass is determined to be 172.5 ± 0.4 (stat.) ± 1.5 (syst.) GeV. 
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1 Introduction 

The top-quark mass is an important parameter of the standard model (SM) of particle physics, 
as it affects predictions of SM observables via radiative corrections. Precise measurements of 
the top-quark mass are critical inputs to global electroweak fits HJI3, which provide constraints 
on the properties of the Higgs boson. 

The top quark constitutes an exception in the quark sector as it decays, primarily to a W boson 
and a b quark, before it can hadronize. Thus, in contrast to all other quarks, the mass of the 
top quark can be measured directly and is currently known with the smallest relative uncer- 
tainty. All measurements of the top-quark mass to date are based on the decay products of tt 
pairs, using final states with zero, one, or two charged leptons. The mass of the top quark has 
been measured very precisely in pp collisions by the Tevatron experiments, and the current 
world average is m t = 173.18 ± 0.56 (stat.) ± 0.75 (syst.) GeV [3|. In the dilepton channel, in 
which each W boson decays into a charged lepton and a neutrino, the top-quark mass has been 
measured to be m t = 170.28 ± 1.95 (stat.) ± 3.13 (syst.) GeV by the CDF Collaboration 0] and 
m t = 174.00 ± 2.36 (stat.) ± 1.44 (syst.) GeV by the DO Collaboration [5]. The combination of 
these two measurements yields a top-quark mass of m t = 171.1 ± 2.1 GeV (3). Measurements 
of m t in pp collisions at \Js = 7 TeV were performed at the Large Hadron Collider (LHC) in the 
dilepton channel by the Compact Muon Solenoid (CMS) Collaboration [6J and in the lepton+jet 
channel, in which one W boson decays into quarks and the other into a charged lepton and a 
neutrino, by the ATLAS and CMS H Collaborations. 

Of all tf decay channels, the dilepton channel has the smallest branching fraction and is ex- 
pected to be the least contaminated by background processes. The dominant background pro- 
cess is Drell-Yan (DY) production. Single top quark production through the tW channel as well 
as diboson production also mimic the dilepton signature but have much lower cross sections. 
The production of multijet events has a large cross section at the LHC, but the contamination 
of the dilepton sample is small as two isolated leptons with high transverse momentum (pj) 
are very rarely produced. The presence of at least two neutrinos in dilepton tt decays gives rise 
to an experimental px imbalance, which allows a further discrimination between background 
and tt events. However, the kinematical system is underconstrained as only the pj imbalance 
can be measured. 

Here we report an update of the measurement of m t performed in dileptonic final states, con- 
taining electrons or muons, with an analytical matrix weighting technique. An alternative 
measurement is performed using a full kinematic analysis. The data samples used in this anal- 
ysis were recorded by the CMS experiment at a centre-of-mass energy of 7 TeV and correspond 
to a total integrated luminosity of 5.0 ± 0.1 fb _1 . 

2 The CMS detector 

The central feature of the CMS detector is a superconducting solenoid, 13 m in length and 6 m 
in diameter, which provides an axial magnetic field of 3.8 T. The bore of the solenoid is outfit- 
ted with various particle detection systems. Charged particle trajectories are measured by the 
silicon pixel and strip subdetectors, covering < <p < 2n in azimuth and \tj\ < 2.5, where the 
pseudorapidity rj is defined as rj = — ln[tan(0/2)], with 9 being the polar angle of the trajectory 
of the particle with respect to the anticlockwise-beam direction. A lead tungstate crystal electro- 
magnetic calorimeter (ECAL) and a brass/ scintillator sampling hadronic calorimeter (HCAL) 
surround the tracking volume; in this analysis the calorimetry provides high-resolution energy 
and direction measurements of electrons and hadronic jets. Muons are measured in drift tubes, 
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3 Simulation of signal and background events 



cathode strip chambers, and resistive plate chambers embedded in the flux-return yoke of the 
solenoid. The detector is nearly hermetic, allowing for pj imbalance measurements in the plane 
transverse to the beam directions. A two-level trigger system selects the most interesting pp 
collision events for use in physics analysis. A detailed description of the CMS detector can be 
found in Ref.0. 

3 Simulation of signal and background events 

The simulation of tt events is performed using the MadGraph [TlOl event generator (v. 5.1.1.0), 
where the generated top-quark pairs are accompanied by up to three additional high-px jets. 
The parton configurations generated by MadGraph are processed with PYTHIA 6.424 |(TT| 
to provide showering of the generated particles. The parton showers are matched using the 
kj-MLM prescription |fl2|. The underlying event is described with the Z2 tune [13] and the 
CTEQ6.6L |fi"4"| set of parton distribution functions (PDFs) are used. The TAUOLA package 
(v. 27.121.5) [15] is used to simulate decays of the x leptons. Events in which the r leptons 
decay to electrons or muons are taken as part of the signal. 

For the reference sample, a top-quark mass of wit = 172.5 GeV is used. Additional samples 
with masses of 161.5 GeV and between 163.5 and 187.5 GeV in steps of 3 GeV are used. Fur- 
thermore, in order to estimate systematic effects in the modelling of dilepton events, simulated 
signal samples using alternative settings of the parameters are also considered. The following 
parameters are varied: the QCD factorisation and renormalisation scale (defined as the squared 
sum of the four-momenta of the primary partons in the event which is transferred dynamically 
in the hard interaction) and the threshold used for the matching of the partons from matrix 
elements to the parton showers. The uncertainty on the choice of the Q 2 or matching scales are 
considered by varying the corresponding nominal value by a factor of two, up and down. 

Electroweak production of single top quarks is simulated using POWHEG (v. 301) fl6l ; Mad- 
Graph is used to simulate W/Z events with up to four jets. Production of WW, WZ, and ZZ 
is simulated with PYTHIA. 

Signal and background processes used in the analysis of tt events are normalised to next-to- 
leading order (NLO) or next-to-next-to-leading order (NNLO) cross section calculations, where 
calculations are available. The production cross section of cr fi = 164+Jgpb computed with 
HATHOR HY7[ [TBI at approximate NNLO is used. The single top quark associated produc- 
tion (tW) cross section is taken to be a m = 15.7 ± 1.2 pb at NNLO 03O. The inclusive NNLO 
cross section of the production of W bosons (multiplied by the leptonic branching fraction of 
the W boson) is estimated to be cr w ^ £v = 31.3 ± 1.6 nb using FEWZ (201 with a Q 2 scale of 
(niw) 2 + E(Px art ° n ) 2 ' wriere m w = 80.4 GeV and p T arton are the transverse momenta of the par- 
tons in the event. The DY production cross section at NNLO is calculated using FEWZ to be 
crz/-y*^u{ m U. > 20 GeV) = 5.00 ± 0.27 nb, where m# is the invariant mass of the two lep- 
tons. In the computation the scales are set using the Z-boson mass niz = 91.2 GeV [21]. The 
normalisation of WW, WZ, and ZZ production is defined using the inclusive cross sections of 
43.0 ± 1.5 pb, 18.8 ± 0.7 pb, and 7.4 ± 0.2 pb respectively (all calculated at NLO with MCFM [22]). 

All generated events are passed through the full simulation of the CMS detector based on 
Geant4 [23 1 . We simulate additional soft Monte Carlo events corresponding to a number of 
collisions distributed as seen in data. 
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4 Event selection 

The tt candidate events are required to contain at least two jets, two energetic isolated leptons 
(electrons or muons), and missing transverse energy (E™ 1SS ) which is defined as the magnitude 
of the pj imbalance vector. Events are selected by dilepton triggers in which two muons, two 
electrons, or one electron and one muon are required to be present. The instantaneous lumi- 
nosity increased significantly during the data taking period thus the lepton pj thresholds were 
increased during the data taking period to keep the trigger rates within the capabilities of the 
data acquisition system. For the dimuon trigger, the pj requirements evolved from 7 GeV for 
each muon to asymmetric requirements of 17 GeV for the highest-px (leading) muon and 8 GeV 
for the second-highest pj muon. For the dielectron trigger, the requirement was asymmetric 
with a threshold applied to the energy of an ECAL cluster projected onto the plane transverse to 
the nominal beam line (Ej). The cluster of the leading electron is required to have Ej > 17 GeV 
and the second-leading electron Ej > 8 GeV. For the electron-muon trigger, the thresholds 
were either Ej > 17 GeV for the electron and pj > 8 GeV for the muon, or Ej > 8 GeV for the 
electron and pj > 17 GeV for the muon. 

All objects are reconstructed using a particle-flow algorithm [24J. The particle-flow algorithm 
combines the information from all subdetectors to identify and reconstruct all particles pro- 
duced in the collision, namely charged hadrons, photons, neutral hadrons, muons, and elec- 
trons. Jets are reconstructed by the anti-fcx jet clustering algorithm |25| with a distance param- 
eter R = 0.5. Jet energy corrections are applied to all the jets in data and simulation [26 j. The 
£miss vec t or i s calculated using all reconstructed particles. 

Events are selected with two isolated, oppositely charged leptons with pj > 20 GeV and 
\t] | < 2.4, and at least two jets with pj > 30 GeV and \tj\ < 2.4. The lepton isolation I re [ is 
defined as the sum of the transverse momenta of stable charged hadrons, neutral hadrons and 
photons in a cone of AR = a/ (A?/) 2 + (A(p) 2 = 0.3 around the lepton track, divided by its trans- 
verse momentum. A lepton candidate is not considered as isolated and is rejected if the value of 
I Te \ is > 0.20 for a muon and > 0.17 for an electron. The two leptons of highest pj are chosen for 
the reconstruction of the top quark candidates. The choice of the jets is different in each anal- 
ysis and is described later. The reconstructed E™ ss of events with same-flavour lepton pairs 
is required to be above 40 GeV to reject DY events. No such selection is applied to e}i events. 
The selected leptons and jets are required to originate from the primary pp interaction vertex, 
identified as the reconstructed vertex with the largest p\ of its associated tracks. Events with 
same-flavour lepton pairs in the dilepton mass window 76 < tnu < 106 GeV are removed to 
suppress the dominant DY production background. Dilepton pairs from heavy-flavour reso- 
nances as well as low-mass DY production are also removed by requiring a minimum invariant 
mass of 20 GeV. A highly efficient b-tagging algorithm based on a likelihood method that com- 
bines information about impact parameter significance, secondary vertex reconstruction, and 
jet kinematic properties, into a b-tagging discriminator, is used to classify the jets [27]. We 
require at least one b-tagged jet in the event. 

The observed number of events is consistent with the expected signal and background yields, 
as shown in Table [l] Simulated events are reweighted to account for differences in trigger, lep- 
ton, and b-tagging selection efficiency between data and simulation. The b-tagging efficiency 
is estimated from a sample of top-quark candidates [28 1 while the probability of tagging light- 
quark jets (mistag rate) is estimated from multijet events ||27|. The lepton selection efficiency 
data-to-simulation scale factors are estimated using dileptons inside the Z-boson mass win- 
dow. The trigger efficiencies are estimated using a data sample collected with a trigger based 
on E™ ss that is weakly correlated with the dilepton triggers and after selecting dilepton events 
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5 Analytical matrix weighting technique 



Table 1: Numbers of observed and expected events in each dilepton channel after all selec- 
tion requirements have been applied. Event yields correspond to an integrated luminosity of 
5.0 fb -1 . The uncertainties quoted correspond to the limited statistics in simulation. The total 
uncertainty associated to the estimates from data of the tf background and DY production are 
included as well. 



Processes 


ee 


epi 


W 


1 b-tagged jet 


tf signal 


598 ± 18 


2359 ± 71 


770 ± 23 


tf background 


10.6 ± 0.3 


101.8 ± 3.1 


15.7 ± 0.5 


Single top 


40.7 ± 1.2 


172.2 ± 5.2 


53.3 ± 1.6 


Drell-Yan 


107 ± 24 


241 ± 27 


143 ± 31 


Dibosons 


11.4 ± 0.3 


39.7 ± 1.2 


13.0 ± 0.4 


Total prediction 


767 ± 30 


2914 ± 76 


995 ± 39 


Data 


817 


2788 


1032 


> 2 b-tagged jets 


tf signal 


1057 ± 32 


4312 ± 129 


1393 ± 42 


tf background 


4.6 ± 0.3 


37.6 ± 1.1 


5.5 ± 0.5 


Single top 


36.8 ± 1.1 


140.6 ± 4.2 


48.2 ±1.4 


Drell-Yan 


38 ± 11 


38.9 ± 4.3 


32 ± 12 


Dibosons 


2.9 ± 0.1 


9.1 ± 0.3 


2.5 ± 0.1 


Total prediction 


1139 ± 34 


4539 ± 130 


1481 ± 43 


Data 


1151 


4365 


1474 



which fulfil the complete event selection criteria. 

The contribution of the DY background is measured using data. For the ee and fip. channels, the 
^out/in method is used |6 1. In this method, the number of DY events counted inside the Z-boson 
mass window in the data is rescaled by the ratio of DY events predicted by the simulation out- 
side and inside the mass window. As contamination from non-DY backgrounds is expected to 
be present in the Z-boson mass window, a subtraction based on data is applied using the eji 
channel scaled according to the event yields in the ee and }iu channels. For the eji channel, 
the DY background yield is estimated after performing a binned maximum-likelihood fit to the 
dilepton invariant mass distribution. The fitting functions are taken from simulation for both 
the signal and background contributions. The contamination from multijet and W+jets back- 
grounds is estimated with a matrix method [29 [, and non-dileptonic tf decays are reweighed 
in the simulation to take these backgrounds into account. This component will be called tf 
background in the following. 

5 Analytical matrix weighting technique 

Since the dilepton channel contains in the final state at least two neutrinos which can not be de- 
tected, the reconstruction of m t from dilepton events involves an underconstrained system. For 
each tf event, the kinematic properties are fully specified by 24 parameters, which are the four- 
momenta of the six particles in the final state: two charged leptons, two neutrinos and two jets. 
Out of the 24 free parameters, 14 are inferred from measurements (the three-momenta of the 
jets and leptons, and the two components of the E™ 1SS ) and 9 are constrained. Two constraints 
arise from demanding that the reconstructed W-boson masses be equal to the world-average 
measured value [21J and one constraint is imposed by assuming the top quark and antiquark 
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masses to be the same f|30l . Further, the masses of the 6 final-state particles are taken as the 
world-average measured values [21 j. This leaves one free parameter that must be constrained 
by using some hypotheses. 

Several methods have been developed for measuring the top-quark mass in the dilepton decay 
channel. We use an improved version of the Matrix Weighting Technique (MWT) [31] that 
was used in the first measurements in this channel |3Tl 1321 . The algorithm is referred to as 
the analytical MWT (AMWT) method. A key improvement with respect to the original MWT 
is the selection of the jets used to reconstruct the top quark candidates. Instead of taking the 
two leading jets (i.e. the jets with the highest pj) , the fraction of correctly assigned jets can be 
increased significantly by using the information provided by b-tagging. Therefore, the leading 
b-tagged jets are used in the reconstruction, even if they are not the leading jets. If there is 
a single b-tagged jet in the event, it is supplemented by the leading untagged jet. The same 
b-tagging algorithm is used as in the event selection. A further improvement is the use of an 
analytical method [33] to determine the momenta of the two neutrinos instead of a numerical 
method. 

In the AMWT, the mass of the top quark is used to fully constrain the tf system. For a given 
top-quark mass hypothesis, the constraints and the measured observables restrict the trans- 
verse momenta of the neutrinos to lie on ellipses in the p x -Py plane. If we assume that the mea- 
sured missing transverse energy is solely due to the neutrinos, the two ellipses constraining 
the transverse momenta of the neutrinos can be obtained, and the intersections of the ellipses 
provide the solutions that fulfill the constraints. With two possible lepton-jet combinations, 
there are up to eight solutions for the neutrino momenta for a given top-quark mass hypoth- 
esis. Nevertheless, in this method, an irreducible singularity that precludes the determination 
of the longitudinal momentum of the neutrinos remains in a limited kinematical region. The 
fraction of events affected by this singularity is below 0.1%, and a numerical method is used to 
determine the solutions in these rare cases [34]. 

The kinematic equations are solved many times per event using a series of top-quark mass 
hypotheses between 100 and 400 GeV in 1 GeV steps. Typically, solutions are found for the neu- 
trino momenta that are consistent with all constraints for large intervals of mass hypotheses. 
In order to determine a preferred mass hypothesis, a weight w is assigned to each solution 11351 : 

w={£f(x 1 )f(x 2 )}p(E* e+ \m t )p(E}-\mt), (1) 

where x, are the Bjorken x values of the initial-state partons, f(x) are the parton distribution 
functions, and the summation is over the possible leading-order initial-state partons (uu, uu, 
dd, dd, and gg). Each term of the form p(E*\m t ) is the probability density of observing a 
massless charged lepton of energy E* in the rest frame of the top quark, for a given m t [35 [: 

Detector resolution effects are accounted for by reconstructing the event 1000 times, each time 
varying the pj, t], and (p of each jet according to the measured detector resolution, and cor- 
recting the E™ 1SS accordingly. For each mass hypothesis, the weights w from all solutions are 
summed. For each event, the top-quark mass hypothesis with the maximum weight is taken as 
the reconstructed top-quark mass w2amwt- Events that have no solutions or that have a max- 
imum weight below a threshold are discarded. This removes 14.6% of the events, and 9934 
events remain in the data, 1550 ee events, 6222 ey. events, and 2110 events. 
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6 Systematic uncertainties 




Generated mass [GeV] Generated mass [GeV] 



Figure 1: Mean mass bias (left) and pull width (right) for different top-quark masses in pseudo- 
experiments for the AMWT method. The red solid line represents the linear fit used to deter- 
mine the correction to apply in order to minimise the residual bias and the blue dashed line 
show the expectation for an unbiased fit. The average pull width for the different top-quark 
masses is 0.99. 

A likelihood C is computed for values of m t between 161.5 and 184.5 GeV, from data in the 
range 100 < Wamwt < 300 GeV. For each value of mt, the likelihood is computed by com- 
paring the reconstructed mass distribution in data with the expectation from simulation. For 
the background, the reconstructed mass distribution of each individual process is added ac- 
cording to its expected relative contribution. Two different templates are used according to the 
b-tag multiplicity of the event, either one b-tagged jet, or two or more b-tagged jets. For the 
DY background, the relative contribution is derived from data in the Z-boson mass window. 
For the other processes, the contributions predicted by the simulation are used. The value that 
maximises the likelihood is calculated after fitting a quadratic function to the — In C values ob- 
tained for all mass points and it is taken as the measurement of m t . Using all the mass points 
in this fit yields pull widths that are closer to unity. 

We determine the bias of this estimate using ensembles of pseudo-experiments based on the 
expected numbers of signal and background events, as shown in Fig.[TJ Given the fit to the data, 
a correction of —0.34 ± 0.20 GeV is applied to the final result to compensate for the residual bias 
introduced by the fit (Fig. [lj left). This correction is obtained from the fit of a linear function 
to the average top-quark masses measured for different mass hypotheses. The width of the 
pull distribution is within 10% of unity for all the mass points, indicating that the statistical 
uncertainties are correctly estimated (Fig. [lj right). 

After correction for the bias, the top-quark mass is measured to be m t = 172.50 ± 0.43 (stat.) GeV. 
The predicted distribution of the reconstructed masses wiamwt for a simulated top quark with 
mass nit = 172.5 GeV, superimposed on the distribution observed in data, is shown in Fig. |2] 
The inset shows the distribution of the — 21n(£/ £ max ) points with the quadratic fit used to 
measure m t . The x 2 probability of the fit is 0.36. 
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CMS 2011 , 5.0 fb 1 at Vs = 7 TeV 




Reconstructed mass [GeV] 

Figure 2: Distribution of the reconstructed mass in data and simulation for a top-quark mass 
hypothesis of 172.5 GeV with the AMWT method. All events used in the analysis are included 
in the distribution. The inset shows — 21n(£/ £ m ax) versus m t with the quadratic fit superim- 
posed. 

6 Systematic uncertainties 

The contributions from the different sources of uncertainty are summarised in Table [2] The un- 
certainty of the overall jet energy scale (JES) is the dominant source of uncertainty on wz t . The 
JES is known with an uncertainty of 1-3%, depending on the pj and rj of the jet 1126 1 . Even in a 
high-pileup regime such as the one observed throughout the 2011 data taking period, the JES 
uncertainty is mostly dominated by the uncertainties on the absolute scale, initial- and final- 
state radiation, and corrections arising from the fragmentation and single-particle response in 
the calorimeter. It has been evaluated for 16 independent sources of systematic uncertainty. 
To estimate the effect of each source on the measurement of m t , the (pi, ?/)-dependent uncer- 
tainty is used to shift concurrently the energy of each jet by ±lc with respect to its nominal 
value, and correcting the E™ 1SS accordingly. For each source, pseudo-experiments are gener- 
ated from simulated event samples for which the JES is varied by the relevant uncertainty, 
and the reconstructed top-quark mass distributions are fitted with the templates derived with 
the nominal JES. The average variation of the top-quark mass is used to estimate the system- 
atic uncertainty. The quadratic sum of the variation for each source is taken as the system- 
atic uncertainty. The uncertainty on pileup corrections to the jet energy calibration (5 sources) 
correspond to a combined uncertainty of 0.53 GeV on wit. Another important contribution is 
the overall data-to-simulation scale calibrated in photon+jet events, yielding a 0.51 GeV un- 
certainty. Other contributions are related to limited knowledge of the single-pion response 
(+Q 3 GeV) and fragmentation models (0.3 GeV) used in the extrapolation as a function of jet pj. 
We also include a time-dependent effect (0.2 GeV) related to variations in calorimeter response 
in the endcaps. Residual eta-dependent corrections based on dijet balance studies (6 sources) 
yield a negligible uncertainty on m t (0.03 GeV). All these sources added in quadrature give a 
combined JES uncertainty of I097 GeV. The final component of JES uncertainty corresponds to 
the uncertainty on the modeling of jet flavour dependence of the jet energy scale (ig 66 GeV) 
which is quoted separately in Table [2] 

The uncertainty due to jet energy resolution is evaluated from pseudo-experiments where the 
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6 Systematic uncertainties 



Table 2: List of systematic uncertainties with their contributions to the top-quark mass mea- 
surement. 



Source 


Am t (GeV) 


Jet energy scale 


+0.90 

— 0.9/ 


b-jet energy scale 


+0.76 
—0.66 


Jet energy resolution 


-l-fl 1 A 




±0 14 


unciustereci 


_i_n 1 9 




+0 05 

-LU.UJ 


Mistag rate 


±0.08 


Fit calibration 


±0.40 


Background normalization 


±0.05 


Matching scale 


±0.19 


Renormalisation and factorisation scale 


±0.55 


Pileup 


±0.11 


PDFs 


±0.09 


Underlying event 


±0.26 


Colour reconnection 


±0.13 


Monte Carlo generator 


±0.04 


Total 


±1.48 



jet energy resolution width in the simulation is modified by ±1<7 with respect to its nominal 
width. The uncertainty on the lepton energy scale is observed to have an almost negligible 
effect on the measurement of m t . The uncertainty in the E™ 1SS scale is propagated to the mea- 
surement of m t after subtracting the clustered (i.e. jet energy) and leptonic components, which 
are varied separately as previously described. This procedure takes into account possible cor- 
relations between the different sources of uncertainty. The scale of the residual unclustered 
energy contribution to the E™ 1SS is varied by ±10% and the corresponding variation of the top- 
quark mass measurement is evaluated from pseudo-experiments. 

The uncertainty due to b-tagging efficiency was evaluated by varying the b-tagging efficiency 
and mistag rates of the algorithm by their respective uncertainties l|27l 1281 . The tagging rate 
was varied according to the flavour of the selected jet as determined from the simulation. This 
affects the multiplicity of b-tagged jets and the choice of the jets used in the reconstruction of 
m t . 

The effect of statistical fluctuations in the templates is estimated by splitting the tt sample 
in four independent subsamples and producing independent templates for each. Pseudo- 
experiments are performed using each new signal template, and the RMS variation of the 
average top-quark mass from each template is taken as an estimate of this uncertainty. The 
uncertainty on the calibration of the fit is added to the systematic uncertainty. The contribu- 
tion from the uncertainty in the ratio between the signal and the background used in the fit is 
evaluated by varying by the corresponding uncertainty the expected number of events. The 
variation of the top-quark mass fit is assigned as a systematic uncertainty. 

The effect due to the scale used to match clustered jets to partons (i.e. jet-parton matching) is 
estimated with dedicated samples generated by varying the nominal matching pj thresholds 
from the default of 20 GeV down to 10 GeV and up to 40 GeV. Effects due to the definition of the 
renormalisation and factorisation scales used in the simulation of the signal are studied with 
dedicated Monte Carlo samples with both scales varied by factors of 2 or \. 
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The uncertainty due to pileup is evaluated from pseudo-experiments where the total inelastic 
cross section used to simulate the pileup is varied within its uncertainty, which is estimated 
to be 8%. The uncertainties related to the parton distribution function (PDF) used to model 
the hard scattering of the proton-proton collisions is evaluated from pseudo-experiments for 
which the distribution of m t was obtained after varying parameters of the PDF by ±lc with 
respect to their nominal values and using the PDF4LHC prescription ffT4l[36ll37| . The differ- 
ences found with respect to the nominal prediction are added in quadrature to obtain the total 
PDF uncertainty. The uncertainties due to the underlying event [13] and the colour reconnec- 
tion [38] are evaluated with dedicated samples. The uncertainties due to the underlying event 
are estimated by comparing two alternative pythia tunes with increased and decreased under- 
lying event activity relative to a central tune. The results for the top-quark mass measured in 
pseudo-experiments using the Perugia 2011 tune are thus compared to the Perugia 2011 mpiHi 
and Perugia 2011 Tevatron tunes [39 [. The difference found between the two samples is taken 
as an estimate of the uncertainty in the modelling of the underlying event in our simulation. 
The Perugia 2011 noCR tune is a variant in which colour reconnection effects are not taken into 
account. The difference in the average top-quark mass, measured with and without colour re- 
connection effects, is taken as the estimate for the colour reconnection systematic uncertainty. 
Finally, the uncertainty due to the modelling of the signal templates by the Monte Carlo genera- 
tor are studied by comparing the results of the pseudo-experiments using the reference sample 
to that from a sample generated with the POWHEG generator. 

7 Measurement with the full kinematic analysis 

An alternative measurement is performed using the KINb method |6l and a tighter event se- 
lection. The jet px is required to be at least 35 GeV and the reconstructed E^ s of ep events is 
required to be at least 30 GeV. These tighter requirements are expected to improve the resolu- 
tion of the method. In KINb, as in the AMWT method, the kinematic equations describing the 
tt system are solved many times per event for each lepton-jet combination. The longitudinal 
momentum of the tt system (p") is used as the extra constraint required to solve the equations. 
The jet px, the E™ 1SS direction, and the p% are varied independently according to their reso- 
lutions in order to scan the kinematic phase space consistent with the tf system. The jet px 
resolution is obtained from the data [26|; the p" description, that is minimally dependent on 
m ir is taken from simulation. The solution with the lowest invariant mass of the tt system is 
accepted if the mass difference between the top quark and antiquark masses is less than 3 GeV. 
The combination of leptons and jets yielding the largest number of solutions is chosen, and the 
mass value 7«KiNb is estimated by means of a Gaussian fit to the distribution of solutions in a 
50 GeV window built around the most probable value. A key point in the method is the choice 
of the jets used to reconstruct the top-quark candidate, favouring jets that have higher value of 
the b-tagging discriminator. Simulations demonstrate that the proportion of events in which 
the jets used for the reconstruction are correctly matched to partons from top quark decays is 
increased significantly with respect to a choice based on the two jets with highest pj. Only 
events with solutions contribute to the m t measurement; in simulation, solutions are found for 
80% of signal events and 70% of background events. 

We use a two-component unbinned maximum-likelihood fit to the wijciNb distribution to miti- 
gate the effect of background and signal events with misreconstructed top-quark masses and 
obtain an estimate of m t . The free parameters of the likelihood are m t and the numbers of sig- 
nal and background events. The main background contribution is from the DY events, which is 
estimated from data using a template fit to the angle between the momenta of the two leptons. 
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8 Summary 



Depending on m t , the signal and background templates may resemble each other; therefore 
the number of background events is constrained by a Gaussian term in the likelihood function. 
The parameters of signal and background templates are taken from simulation and fixed in the 
fit. The signal shape is obtained with a simultaneous fit of simulated tt samples to a Gaussian 
plus Landau function template with parameters that are linear functions of m t . Separate tem- 
plates are used for the four samples corresponding to the same or different flavour dileptons 
with one or two and more b-tagged jets. In each category the backgrounds are added in the 
expected proportions. The expected distribution from DY events is determined from data near 
the Z peak (76 < Mm < 106 GeV) for same-flavour dileptons. From simulation, the template 
obtained near the Z peak is expected to describe well DY events in the signal region. In the 
case of different-flavour dileptons we estimate the contribution from DY events using a data 
sample of Z — > jiji, by replacing the muons with fully simulated decays of t leptons [40] and 
applying the event selection and top-quark mass reconstruction. For single top quark, diboson, 
and other residual backgrounds the templates are taken from simulation. 

The fit is performed separately for same- and different-lepton flavour events with either one or 
at least two b-tagged jets using an unbinned likelihood method, where the inputs are the mass 
value returned by the KINb method in the data, and the probability density function for signal 
and background. The data in the range 100 < mjciNb < 300 GeV is used in the fit. Figure [3] 
(inset) shows the variation of — 21n(£/£ max ) as a function of m t , for the different categories 
individually and for all categories combined. For each event category the corresponding like- 
lihood is maximised, yielding an estimate of the top-quark mass value as well as the expected 
numbers of signal and background events. The result of the fit for the category of events with 
the smallest background contamination (eji events with at least two b-tagged jets) is shown in 
Figure |3j 

The expected contamination from background events and the result obtained from the fit in 
each category agree well. A combined unbinned likelihood is constructed in order to extract 
the final measurement of m t from data. To minimise any residual bias resulting from the pa- 
rameterisations of the signal and background mjQNt, distributions, pseudo-experiments are per- 
formed using simulated dilepton events generated with different m t values. The resulting m t 
distributions are used to calibrate the parametrisation of the signal template. We find an aver- 
age bias on m t of 0.4 ± 0.2 GeV, which we use to correct our final value. We assign the envelope 
of the residual bias (0.2 GeV) as the systematic uncertainty associated with the fit. 

Other sources of systematic uncertainty are similar and fully correlated with those in the AMWT 
analysis. We observe however that the KINb method is affected by larger uncertainties com- 
pared to the AMWT method, reflecting the fact that the mass resolution is slightly poorer. 
The degradation of the resolution is related to the fact that a choice is made for the lepton- 
jet assignment in the event and that there is no reweighting of the solutions found based 
on any expectation for the kinematic properties, such as polarization effects which are in- 
trinsically modelled by Eq. [2] We find no improvement in combining the AMWT and KINb 
given the difference in statistical uncertainty achieved and the dominance of the correlated 
systematic uncertainties. The KINb analysis is thus used as a cross-check, and we measure 
m t = 171.8 ± 0.6 (stat.) ± 2.2 (syst.) GeV, in agreement with the AMWT measurement. 

8 Summary 

In summary, a measurement of the top-quark mass from tt decays to dilepton final states is 
presented, using a data sample corresponding to an integrated luminosity of 5.0 fb _1 recorded 



11 



> 450 r 




CMS 201 1 , 5.0 fb 1 at = 7 TeV 



CD 



2 4oo ;_ 



C 



9? 300^ 



LU 



250^ 



200- 



150 



100 



50- 



nt r-T 

TOO 120 140 160 180 200 220 240 260 280 300 

Reconstructed Mass [GeV] 



Figure 3: Result of the fit to eji events with at least two b-tagged jets using the KINb method. 
The inset shows the variation of the likelihood used to extract the top-quark mass, in the dif- 
ferent event categories and for all channels combined. 

by the CMS experiment at y^s = 7 TeV. The measurement yields m t = 172.5 ± 0.4(stat.) ± 
1.5 (syst.) GeV. An alternative measurement gives a consistent result. With respect to the previ- 
ous measurement in the dilepton channel performed by CMS on the 36 pb -1 data collected in 
2010 [6j, the systematic uncertainty could be reduced substantially by improved understanding 
of the effect of pileup, underlying event and the uncertainty on the JES. To date, this measure- 
ment is the most precise determination of the top-quark mass in the dilepton channel. 
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